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Abstract. The Einstein Equivalence Principle (EEP) is one of the foundations of 
the theory of General Relativity and several alternative theories of gravitation predict 
violations of the EEP. Experimental constraints on this fundamental principle of nature 
are therefore of paramount importance. The EEP can be split in three sub-principles: 
the Universality of Free Fall (UFF), the Local Lorentz Invariance (LLl) and the Local 
Position Invariance (LPl). In this paper we propose to use stable clocks in eccentric 
orbits to perform a test of the gravitational redshift, a consequence of the LPL The 
best test to date was performed with the Gravity Probe A (GP-A) experiment in 
1976 with an uncertainty of 1.4 x lO”"*. Our proposal considers the opportunity of 
using Galileo satellites 5 and 6 to improve on the GP-A test uncertainty. We show 
that considering realistic noise and systematic effects, and thanks to a highly eccentric 
orbit, it is possible to improve on the GP-A limit to an uncertainty around (3—4) x 10“® 
after one year of integration of Galileo 5 and 6 data. 
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1. Introduction 

The classical theory of General Relativity (GR) is the current paradigm to describe 
the gravitational interaction. Since its creation in 1915, GR has been conhrmed 
by experimental observations. Although very successful so far, it is nowadays 
commonly admitted that GR is not the ultimate theory of gravitation. Attempts 
to develop a quantum theory of gravitation or to unify gravitation with the others 
fundamental interactions lead to deviations from GR. Moreover, observations requiring 
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the introduction of Dark Matter and Dark Energy are sometimes interpreted as a hint 
that gravitation presents some deviations from GR at large scales. 

GR is built upon two fundamental principles. The hrst principle is the Einstein 
Equivalence Principle (EEP) which gives to gravitation a geometric nature. More 
precisely, the EEP implies that gravitation can be identihed with space-time curvature 
which is mathematically described by a space-time metric If the EEP postulates 
the existence of a metric, the second principle of GR specihes the form of this metric. 
In GR, the metric tensor is determined by solving the Einstein held equations which 
can be derived from the Einstein-Hilbert action. 

From a phenomenological point of view, three aspects of the EEP can be tested [T] : 
(i) the Universality of Free Fall (UFF), (ii) the Local Lorentz Invariance (LLI) and (iii) 
the Local Position Invariance (LPI). The UFF has been tested experimentally at the 
very impressive level of 10“^^ with torsion balances (see e.g. 12112 ]) and with Lunar Laser 
Ranging measurements [1] . One way to test the LLI is to search for anisotropies in the 
speed of electromagnetic interactions, usually parametrized with the c^-formalism [1]. 
It is interesting to mention that GPS data have provided a good test of the LLI j5], even 
if a better constraint is now given by Ives-Stilwell type experiments [B]. The last part 
of the EEP, the LPI is tested by constraining space-time variations of the constants of 
Nature (see e.g. [U [HI [HI HD] ) or by redshift tests [HI [121 [13] • 

In addition to the EEP tests, the second principle of GR (the form of the metric) is 
thoroughly tested by different observations in the Solar System: deflection of light jH] , 
planetary ephemerides HSlttniE], radioscience tracking data of spacecraft HElEHlEn], 
etc (for a larger review, see [U [2T]L 

Tests of the EEP are of prime importance since violations thereof are predicted by 
attempts to develop a quantum theory of gravity or a theory unifying all fundamental 
interactions and in several models of Dark Energy (see e.g. [22] and references therein for 
some examples). The level of the EEP violation depends highly on the theoretical model 
considered (see for example [23]) while the different types of EEP tests (universality of 
free fall, redshift test, constancy of the constants of Nature) are sensitive to different 
fundamental parameters 1211 and are very complementary (see also the discussion in 
section V.A. of [23]L 

As mentioned above, a gravitational redshift experiment tests the LPI [1]. The 
most precise test of the gravitational redshift to date has been realized with the 
Vessot-Levine rocket experiment in 1976, also named the Gravity Probe A (GP-A) 
experiment [13]. The frequency differences between a space-borne hydrogen 

maser clock and ground hydrogen masers were measured thanks to a continuous two-way 
microwave link. The gravitational redshift was verihed to 1.4 x 10“^ accuracy [13]. The 
future Atomic Glock Ensemble in Space (AGES) experiment, an ESA/GNES mission, 
planned to fly on the ISS in 2017, will test the gravitational redshift to around 2—3 x 10“® 
accuracy [23]. Furthermore, other projects like STE-QUEST [25] propose to test the 
gravitational redshift at the level of 10“^, and observations with the RadioAstron 
telescope may reach an accuracy of the order of 10 -'■ 122 ]. Finally, it has been previously 
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suggested in [2^ to use Galileo satellites for such a test. 

In this paper, we will show how clocks on-board GNSS satellites can be used to 
perform an improved test of the gravitational redshift if they are placed on an elliptical 
orbit around Earth. As an application, we will consider the case of the satellites Galileo 
5 and 6. These satellites were launched on August, 30th 2014 and because of a technical 
problem, the launcher brought them on a wrong, elliptic orbit. An elliptic orbit induces 
a periodic modulation of the gravitational redshift [22] while the good stability of recent 
GNSS clocks allows to test this periodic modulation to a very good level of accuracy. 
The Galileo 5 and 6 satellites, with their large eccentricity and on-board H-maser clocks, 
are hence perfect candidates to perform this test. Gontrary to the GP-A experiment, it 
is possible to integrate the signal on a long duration, therefore improving the statistics. 
In this paper, we assess the sensitivity of the gravitational redshift test that can be 
performed with Galileo satellites 5 and 6. 

In Sec. we describe the experimental data needed in order to perform the test 
with Galileo satellites. Then, we introduce in Sec. the simulation of the satellites 
orbits, clocks and signal. In Sec. we estimate the test sensitivity with two different 
statistical methods to estimate the averaging of random noise in the signal. Finally, in 
Sec. we present a detailed study of the different systematics that need to be modeled 
in our analysis and we show how to decorrelate them from the gravitational redshift 
signature. 

2. Experimental data 

The experiment proposed here requires an accurate knowledge of the frequency of the 
satellite clock as it orbits the Earth. These data are made available by several analysis 
centers (AGs) of the International GNSS Service (IGS, [30]) in the framework of the 
Multi-GNSS-EXperiment (MGEX, [3l]). MGEX applies the standard IGS processing 
chain to new GNSS (Galileo, Beidu) and regional services (e.g., QZSS for Japan). 
All active Galileo satellites orbits, clock corrections, and inter-system biases are then 
available for public use through the MGEX product directory at the IGS data center 
GDDIS, which is updated on a regular basis. 

Since each AG is following a different strategy, as an example we shall focus on 
the processing performed at the GODE (Genter for Orbit Determination in Europe) 
analysis center [32]. In this particular case, the whole solution is based on a two steps 
process. First, satellite orbits are determined using a double difference network solution. 
This technique allows to eliminate or reduce several biases by combining observations 
among a pair of receivers and a pair of satellites. All clock errors are eliminated in the 
combination while the remaining parameters (including ephemerides, atmospheric and 
antenna parameters, etc...) are estimated for all stations and satellites involved. The 
latter are then used as a framework for the clock solution, which is based on a zero- 
difference processing. For more general information about the IGS and its products and 
data centers, the reader can refer to [3U] . 
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Figure 1. (left) MDEV for Galileo In-Orbit Validation (lOV) (GALll-12-19-20) and 
for GPS Block II-F (PRN GPSOl-25) from [34]. Superimposed in blue is the MDEV 
of the simulated clock noise; (right) Gravitational redshift signal with a = 0. 


For the test proposed in this paper, the orbit solution of Galileo 5 and 6 satellites 
would be used to calculate the behaviour of the on-board clocks and the gravitational 
redshift as predicted by GR |2I]. The latter would then be compared to the clock 
solution from the IGS processing to recover any violation of the LPI, as detailed in 
Sec. 1^ Several conditions need to be satished by the data to assure the success of the 
experiment. The availability of continuous orbit and clock parameters is paramount and 
is assured by the IGS guidelines for the clock and orbit products. Orbits and station 
coordinates should be unaffected by an eventual violation of the gravitational redshift. 
This is assured as a result of the processing outlined above, since both are derived by a 
double-difference solution which eliminates all clock parameters. Finally, the reference 
clock for the solution is always chosen among ground station H-maser clocks [33], so 
that it is unaffected by an eventual LPI violation. 

3. Simulation: orbits, clocks and signal 

Orbital parameters of Galileo satellites 5 and 6 used for the simulation of the present 
work are summarized in table [Tj Their orbits have been recently circularized, from 
eccentricity e = 0.2330 to e = 0.1561. However, we will consider both the initial and 
the hnal orbits, in order to compare the sensitivity of the redshift test to different 


Name 

e 

a (km) 


Galileo 5&6 Final Orbit 

0.1561 

27977 

49.7212 

Galileo 5&6 Initial Orbit 

0.2330 

26192 

49.7740 


Table 1. Orbital parameters of satellites Galileo 5 and 6 [3S|. Their orbits were 
circularized so that they can be used for positioning. Both orbits are now very similar. 
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eccentricities. Orbits are calculated by solving the Kepler equation for a duration of 
two years, thus ensuring the numerical stability of the solution with time. 

We use the clock solutions provided by CODE to simulate a realistic random noise 
for Galileo satellite clocks. Fig. (left) shows the Modified Allan Deviations (MDEVs) 
for several GPS and Galileo clocks |31] around Day Of Year (DOY) 100/2013, as well 
as the MDEV of the simulated clock noise (in blue). For an integration time of 1000 s, 
the MDEVs are between 3 x 10“^"^ and 7 x 10“^^, and the MDEV of the simulated clock 
noise is around 5 x 10“^^. The MDEVs decrease like which is a characteristic 

of white frequency noise. There is a bump at about 21000 s due to systematic effects, 
which will be analysed in Sec. Moreover, we add a flicker noise to the simulated clock 
noise at 8 x 10“^® in the MDEV (flat part of the blue curve on Fig. left). Although 
Fig. 0 (left) shows no sign of an increase of noise at large averaging times (indicative 
of random walk frequency noise or frequency drift), on ground measurements of Galileo 
H-maser clocks have shown such behaviour [36]. The latter is anyway not exceeding 
5 X 10“^^ in MDEV for measurements up to 10^ s averaging time. At the frequency of 
our expected periodic signal (2.1 x 10“® Hz), the corresponding noise is about an order 
of magnitude below our noise model, and we will therefore neglect it in our analysis. 

The ideal signal y is the gravitational part of the relative frequency difference 
between a ground clock g and the satellite clock s: 


y = 


GM ( 1 


( 1 ) 


where GM is the standard gravitational parameter of Earth, c is the speed of light in 
vacuum, and Vg and are respectively the norm of the position vectors of the ground 
clock and of the on-board clock. We neglect the noise coming from the ground clock, 
and we are interested only in the variable part of y. Therefore one can remove the first 
term in ([^, which is constant. 

We use a simple phenomenological model for a possible violation of the gravitational 
redshift, characterised by a parameter a which is zero in general relativity ID EH: 


^( a ) =-(1 + a )^ . ( 2 ) 

The GP-A limit corresponds [13] to cTq, = 1.4 x 10“^, where cTq, is the uncertainty on the 
estimation of the parameter a. 

The simulated signal y is shown in Fig. (right). Its amplitude is around twice 
larger for Galileo satellites 5&6 initial orbits than for Galileo 5&6 final orbits. The 
asymmetry between the upper part and the bottom part of the curves is a characteristic 
of the elliptical orbit. 


4. Estimation of the statistical sensitivity of the test 

We study the sensitivity of the gravitational redshift test performed with Galileo 5 and 6 
through two very different methods: the Fast Fourier Transform (EFT) method and the 
Linear Least-Square (LSQ) method. 




Test of the gravitational redshift, application to Galileo satellites 5&6 


6 


a 

.2 

'o 

■> 

I—I 

Ph 

1-1 

n 

o 


OJ 

o 



Figure 2. Statistical sensitivity of the gravitational redshift test with respect to the 
duration of the experiment. We show the results for the two methods FFT and LSQ, 
and for the initial and final orbits of Galileo satellites 5&6. In green is the accuracy of 
the Gravity Probe A experiment for reference. 


The sensitivity of the gravitational redshift test can be calcnlated as the inverse of 
the Signal-to-Noise Ratio (SNR) p. The SNR can be maximized with matched hltering, 
i.e. by nsing a hlter optimized for the searched signal and the considered noise. Then 
the SNR is given by [37|: 


P 


2 



Seif) 


(3) 


where y and e are the gravitational redshift signal and the clock noise, respectively, J^[y] 
is the Fonrier transform of y, and is the Power Spectral Density (PSD) of the clock 
noise. We implemented a nnmerical version of this method with the following steps: 
hrst we generate 100 different seqnences of clock noise; we calcnlate the PSD of each 
seqnence nsing a fast Fonrier transform algorithm and average the 100 different PSD; 
hnally we compnte the Signal-to-Noise Ratio (SNR) ([^ with the mean PSD of the clock 
noise, by integrating in the freqnency domain [—Fs/2,Fs/2], where Fg is the sampling 
rate of the signal. This is done for several different dnrations of the experiment. 

The second method implemented consists in using a linear least-square (LSQ) £t 
combined with a Monte-Carlo method, where we find the minimum of the merit function 

N 

^=1 


( 4 ) 
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with respect to a and A. In Eq. Q is the nnmber of simnlated observables, A 
is a constant to be estimated and (ej)j=i.jv is one seqnence of simnlated clock noise, 
while (y(ti) + €i) corresponds to the simulated observable at time U. The constant A is 
introduced to remove a constant frequency bias which cannot be measured accurately. 
Instead of fitting A, it is also possible to remove the mean of y{a) from the model, so 
that the fitted model has a zero mean, like the noise e. We implemented a numerical 
version of this method, with the following steps: we generate 100 different sequences of 
clock noise (we take the same ones as for the EFT method); we estimate a for each clock 
noise sequence with the linear LSQ method, for different durations of the experiment; 
finally, for each duration, we calculate the mean of the 100 obtained values of a and 
their standard deviation (Jq. The standard deviation corresponds approximately to 68% 
of the obtained values around the mean. 

Using [38], it can be shown that for high SNR and in the case of a one parameter 
determination, the inverse of the SNR found with matched filtering method is equal to 
the standard deviation of a found with the LSQ method, i.e. cxa = p~^. 

Fig-i shows the statistical sensitivity of the gravitational redshift test, aa, with 
respect to the duration of the experiment for the two different methods presented above 
and for the final and initial orbits of Galileo satellites 5 & 6. Also, we verify that the 
two methods give the same estimate, thus giving us strong confidence in the reliability 
of our analysis. This result is very promising, as it can be seen that even with the less 
eccentric orbit, the GP-A limit could be attained in less than a month. However, this 
result has to be tempered with a study of systematic effects. 

5. Systematic effects 

All systematic effects (i.e. errors that are due to unmodeled deterministic effects) that 
mimic the gravitational redshift signal may induce a bias in the estimation of a, i.e. 
a “fake” violation of the gravitational redshift. This will ultimately limit the useful 
duration of noise averaging presented in Fig. The bias on the estimated value of a 
resulting from unmodeled systematic effects will depend on their magnitude, frequency 
and phase. We classify systematic effects in four classes: (i) effects acting on the 
frequency of the reference ground clock; (ii) effects on the links, e.g., mismodeling 
of ionospheric or tropospheric delays, variations of receiver/antenna delays, multipath 
effects, etc... ; (iii) effects acting directly on the frequency of the space clock, e.g., 
temperature and/or magnetic field variations on board the Galileo satellites ; (iv) orbit 
modelling errors. The latter directly translate into variations of the measured clock 
frequencies because the clock solutions are obtained from the one-way signals (so called 
“zero differences”), hence an orbit error of Ax along the line of sight leads directly to 
an error in the clock phase estimation of ~ Ax/c. 

Effects on the ground reference clock are likely to be negligible with respect to the 
others, and are furthermore likely to be only weakly correlated with the signal we are 
searching for shown in Fig. (right). They will therefore be neglected in this analysis. 
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Effects on the link could be significant for certain stations (e.g. under unfavourable 
multipath environment or atmospheric conditions). However they are expected to 
vary with elevation and azimuth, and local meteorological conditions, none of which 
are correlated with the expected signal. Note also that the main signal frequency is 
2.1 X 10“® Hz whereas the 2nd harmonic of any diurnal effect (temperature, humidity 
etc...) is 2.3 X 10“® Hz, which indicates that such effect will decorrelate from the signal 
with as little as 6 days of data. Finally, some more significant averaging can be expected 
given the large number of observing stations. 

The environmental sensitivity of the Galileo passive hydrogen masers (PHM) has 
been characterized on ground (see e.g., [SB])- Temperature sensitivity is < 2 x 
and magnetic field sensitivity < 3 x 10“^^/G. Both Galileo 5 and 6 are equipped with 
on-board temperature sensors, but no in-situ temperature measurements are publicly 
available. From Fig. (right) we see that the expected peak to peak variation of 
the redshift signal is about 5 x 10“^^. Given the expected statistical sensitivity of 
our test (cf. Fig. we want any systematic effect to remain 2 x 10“® below that 
level, i.e. < 10“^^. This implies that temperature variations at the clock that are 
correlated with the signal need to be < 0.05 K peak to peak. Whether this is achieved 
in orbit remains to be seen once temperature measurements are made available. As an 
example, [22] mention that for the GPS block HF “the rubidium clocks of the Block 
HF satellites are mounted on a thermally isolated sub-panel (along with the frequency 
distribution unit) for which a peak-to-peak temperature variation of less than 0.5 K has 
been observed in SVN62 on-board measurements. In addition, the new clock benefits 
from a baseplate temperature controller (BTC), which further reduces the sensitivity to 
external temperature variations by up to a factor of 50.” This indicates that peak to 
peak variations < 0.05 K are realistic. Additionally, significant decorrelation between 
temperature effects and our signal can be expected, as the former depend on Sun 
exposure of the satellite and thus have annual spectral components, which are absent 
in the signal. 

To evaluate the influence of the magnetic effect on our determination of a we 
model the Earth magnetic field as a dipole with the field amplitude at the location of 
the satellite given by 


B 


Bo 


Vl + 3 cos^ 6 , 


(5) 


where Bq = 0.312 G, Re is the mean radius of the Earth and 6 is the angle between 
the position vector of the satellite and the geomagnetic North pole. We then use the 
clock sensitivity coefficient to add a systematic effect esys(t) to our signal, where the 
time dependence is given by ([^ and by the satellite orbit. We determine a from the 
simulated signal Y using the LSQ method (sec. |^, such that 


^ (f) — y(t) + ^sjsit) + ^(t), ( 6 ) 

where e is the random clock noise. Fig. shows the obtained value of a (the bias on 
a induced by the systematic effect) as a function of observation time together with the 
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Figure 3. Statistical uncertainty (blue circles) on the LPI violation and bias induced 
by uncorrected magnetic effect (black line) as a function of observation time, for 
Galileo 5&6 final orbits. 


obtained uncertainty from the LSQ-Monte Carlo method. We see that for observation 
times up to 450 days the magnetic effect remains below the statistical uncertainty. For 
longer integration times it will become limiting at ~ 2.7 x 10“®. This is a worst case 
estimation as it is of course possible to correct at least partly for the magnetic effect 
and/or simultaneously fit the magnetic effect and a (c.f. the discussion of the solar 
radiation pressure effect below) leading to a better result. 

Finally, we look at effects related to orbit modelling errors. The radial part of 
the resulting orbit error contributes directly to the clock estimation error. Indeed, it 
has been shown that the clock offset is correlated with satellite laser ranging (SLR) 
residuals, which can be considered as a measure of the radial orbit error [10]. The 
clock error also shows a dependency with the Sun elevation angle /3, which is the angle 
between the satellite orbital plane and the direction of the Sun. This clearly indicates 
systematic effects which are linked to the direction of the Sun, as for a mismodeling of 
Solar Radiation Pressure (SRP) [lO]. Indeed, the use of an enhanced SRP model for 
orbit determination can divide by a factor of four the effect of SRP mismodeling on the 
clock estimation. The resulting systematic effect on the clock estimation is evidenced 
as a bump at half-orbital time scale in the Modihed Allan Deviation (MDEV) of Fig. 
(left). 

Consequently, we will consider here a class of systematic effects which shows a 
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Figure 4. Statistical uncertainty of the determination of the gravitational redshift 
in the presence of solar radiation pressure, as a function of observation time, for 
Galileo 5&6 final orbits, and for six different occurrences of the simulated clock noise. 


dependency with the Sun direction, as well as with the mean anomaly of the satellite. 
For the SRP effect in particular we use a model that reproduces the annual amplitude 
variation (cf. Fig. 3 of [10]) and has a frequency of n + cua, where n is the mean motion 
of the satellite and LVa = 27r/year: 

^sys(i) = ^(1 + B(cos(uat + 0i) — 1)) sin((n + Ua)t + 02) , (7) 

with A, B,(j)i,(j )2 constant. We simulate data with A = 5 x 10“^^, R = 0.45, 0i = 0 
that reproduce the results of [ID], and different values of 02- We then use a Bayesian 
approaclijl] with a slice sampling method [H] assuming normal distributions and flat 
priors to simultaneously determine the four parameters of the SRP model ([^ and the 
redshift parameter a, as well as their respective uncertainties. The result is shown in 
Fig. 1^ for six different occurrences of the simulated clock noise. We hnd that for data 
spanning more than 1 year, and a “worst case” choice of 02 the resulting value of a is 
unbiased with a statistical uncertainty of ~ 4 x 10“®, slightly increased with respect 
to the uncertainty when determining a with no SRP effect present in the observable 
(c.f. Fig.g). After 200 days of integration time, the six occurrences of noise converge 
to the same values which indicates that the uncertainty is limited by the residual SRP 
systematic effect. In our estimates we have assumed the “standard” SRP corrections 


f See for example 

http://fr.mathworks.com/help/stats/examples/ 

bayesian-analysis-for-a-logistic-regression-model.html 
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rather than the improved ones discussed in mi. The latter should allow a reduction 
of the effect by a factor 2-4 [40], so our estimates are likely to be again only an upper 
limit. 


6. Conclusion 

In this communication, we have shown how stable clocks placed in eccentric orbits can 
provide a powerful test of the gravitational redshift. These tests are of high scientihc 
relevance, as many alternative theories of gravitation predict a EEP violation at some 
level of accuracy (see for examples |l2l [23l |22l 113|). In particular, we have shown that 
the Galileo 5 and 6 GNSS satellites can improve on the GP-A (1976) limit on the 
gravitational redshift test, down to an accuracy of (3 — 4) x 10“® with at least one year 
of data. 

Following our study, we recommend the data of Galileo satellites 5 and 6 to be 
made available for at least one year, in order to be able to decorrelate systematic 
effects from the gravitational redshift signal and to reach sufficiently low statistical 
uncertainties. Some data is already being collected by the IGS but ideally it should 
be complemented by housekeeping data such as on-board temperature, which would 
provide further conhdence in the estimation of systematic effects. 
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